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PREFACE 


This  final  technical  report  was  authored  by  Dr.  T.  Vessel  Berg  and  pre- 
pared by  Electron  Dynamics  Division,  Hughes  Aircraft  Company,  Torrance, 
California  on  Contract  F30602-77-C-0221,  for  Rome  Air  Development  Center, 
Grlfflss  Air  Force  Base,  New  York.  It  summarizes  the  results  of  a large 
signal  theory  for  traveling  wave  tubes,  based  on  polorlzatlon  variables, 
that  forms  the  basis  for  a new  large  signal  computer  program  for  TWT. 

This  effort  has  been  Initiated  In  October  1977  and  was  completed  In 
November  1978.  Henry  Friedman  (OCTP)  was  RADC  Project  Engineer. 
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EVALUATION 


A majority  of  TWT  designs  are  developed  through  an  iteration  of 
design  and  test  techniques.  This  results  in  a significant  time  delay 
from  concept  to  the  practice  and  larger  development  costs  due  to  the 
"cut  and  try"  techniques  used.  In  an  effort  to  alleviate  this  problem, 
relatively  sophisticated  ^nd  Company  Proprietary)  computer  codes  were 
generated  to  reduce  turn  around  time  and  design  costs.  Since  most  of 
these  codes  were  based  on  Lagrangian  and  Eulerian  formulations,  simula- 
tions were  quite  lengthy  and  costly.  This  report  describes  a more  cost 
effective  technique  which  gives  good  agreement  with  experimental  TWT 
data.  Significant  computation  time  savings  have  resulted  from  the  use 
of  the  polarization  model  described  herein,  allowing  the  designer  to 
practically  interact  with  the  resultant  simulation.  Further,  by  the 
availability  of  this  code  to  Government  agencies,  it  is  hoped  that  this 
code  will  enable  more  quantitative  evaluation  of  future  TWT  designs. 


Project  Engineer 


1.0  INTRODUCTION 


This  report  presents  the  basis  of  a large-signal  theory  for  traveling 
wave  tubes  assuming  a non-relativistlc  one-dimensional  beam.  The 
theory  is  expressed  in  terms  of  polarization  variables  which  have  not 
been  used  in  large-signal  analyses  of  TWT's  before,  but  can  be  shown  to 
Introduce  simplifications  in  describing  certain  non-linear  effects  in 
the  beam  dynamics  [1,2]. 

The  basic  idea  of  the  present  formulation  can  perhaps  best  be  explained 
by  comparison  with  Lagranglan  and  Eulerian  type  formulations  which  are 
well  known  in  large-signal  TVTT  calculations.  These  latter  formulations 
are  characterized  by  non-linear  beam  dynamics  but  linear  coupling 
between  the  beam  and  the  circuit  field. 

On  the  other  hand  the  polarization  description  is  characterized  by  the 
unique  feature  that  the  dynamic  equations  are  linear  even  in  the  large- 
signal  domain,  whereas  the  coupling  to  the  circuit  field  is  non-linear. 
Hence,  in  the  polarization  model  the  non-linear  problems  are  absent  in 
the  beam  dynamics,  but  they  reappear  in  the  coupling  to  the  circuit 
field,  and  to  a lesser  extent,  in  the  space  charge  field. 

The  basic  dynamic  rf  variables  are  the  displacement,  or  polarization, 
and  the  velocity  at  the  displaced  position.  It  can  be  shown  that  the 
harmonics  of  these  variables,  for  a specified  modulation  in  the  beam, 
decay  faster  than  the  harmonics  of  the  rf  current  and  rf  space  charge 
used  in  the  common  large  signal  descriptions  [3].  This  feature,  which 
is  unique  for  the  polarization  model,  has  the  Important  practical 
implication  that  a smaller  number  of  beam  harmonics  is  sufficient  to 
describe  the  large  signal  beam  dynamics.  In  the  present  work  the 
highest  harmonic  number  Included  is  four,  although  the  basic  formula- 
tion is  for  an  infinite  number  of  harmonics. 


2.0  LARGE  SIGNAL  THEORY  WITH  POLARIZATION  VARIABLES 


In  this  chapter  the  basic  forms  of  the  circuit  equations  and  electronic 
equations,  expressed  In  polarization  variables,  are  presented.  Let  us 
first  consider  the  circuit  equations. 

2.1  CIRCUIT  EQUATIONS 

The  equivalent  circuit  used  in  the  analysis  is  shown  in  Figure  2-1. 

The  inductance  L,  capacitance  C,  and  resistance  R,  all  measured  per  unit 
length,  are  specified  as  a function  of  frequency,  so  that  each  harmonic 
frequency  component  has  its  own  equivalent  circuit.  The  circuit  param- 
eters are  related  to  the  actual  circuit  parameters  by  the  relations: 


L(u))C(a))  = — 

Vc(a)) 


(2-1) 


L((>)) 

C(a)) 


(2-2) 


where  Vj,((o)  is  the  phase  velocity  and  Z^(a>)  the  impedance  of  the  circuit, 
both  referred  to  a particular  frequency  to.  Note  that  the  impedance 
Zc(t^)  is  specified  as  the  impedance  at  the  circuit  radius  and  not  at  the 
location  of  the  beam. 


With  these  definitions  the  circuit  equations  take  the  form: 


^ _1 ^ 

3z  3t 


31 


in 


dz 


(2-3) 


3V  C 31 


(2-4) 


2 


Figure  2- 


Equivalent  circuit  diagram  of  the 
circuit  interacting  with  the  beam. 
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where  I and  V are  the  circuit  current  and  circuit  /oltage,  respectively. 

The  induced  current  into  the  circuit  from  the  beam  is  I.  . 

in 

Two  important  points  should  be  emphasized  in  connection  with  the 
circuit  equations  (2-3)  and  (2-4).  The  first  point  has  to  do  with  the 
choice  of  rf  variables  of  the  circuit,  which  are  the  circuit  current  I 
and  the  circuit  voltage  V.  The  choice  of  voltage  V rather  than  elec- 
tric field  E (which  is  often  used  in  TWT  analyses)  is  dictated  from 
the  requirement  that  the  rf  variables  should  be  continuous  across 
intersections  between  two  adjacent  circuits  characterized  by  different 
circuit  parameters.  Such  configurations  exist  in  the  usual  TWT  con- 
sisting of  the  input  section,  the  attenuator  section,  and  the  output 
section.  It  is  easy  to  demonstrate  that  the  electric  field  E is  discon- 
tinuous across  such  intersections,  and  therefore  is  unacceptable  as  rf 
variable.  Since  E =■  -3V/3z,  and  I is  continuous,  it  follows  directly 
from  (2-4)  that  E is  discontinuous  across  an  intersection  characterized 
by  unequal  circuit  parameters  Z^,  and  R on  the  two  sides  of  the 
intersection.  Since  continuity  of  variables  is  essential  in  the  match- 
ing procedure,  it  follows  that  the  correct  choice  of  rf  circuit  varia- 
bles is  I and  V,  rather  than  I and  E. 

The  second  point  to  be  noticed  in  connection  with  (2-3)  is  that  the 

Induced  current  I.  in  the  circuit  from  the  beam  is  evaluated  somewhat 
in 

differently  from  the  usual  way. 

2.1.1  The  Driving  Current  in  the  Circuit  Equation 

In  the  usual  theory  of  TWT's  the  driving  current  is  taken  to  be  the 
beam  current  itself,  whereas  the  circuit  impedance  is  reduced  by  an 
appropriate  factor,  thereby  essentially  placing  the  circuit  at  the 
position  of  the  electron  beam.  In  the  present  theory  a different  and 
more  consistent  procedure  is  followed,  in  which  the  finite  spacing 
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between  circuit  and  beam  Is  accounted  for  by  a coupling  coefficient, 
much  In  the  same  way  as  usually  done  In  klystrons. 

The  Induced  current  flowing  In  the  circuit  Is  obtained  from  the 
Integral: 


KZo*  ~ f ^b^*’ 

X 

In  this  equation,  t)  Is  the  beam  current  density,  f(x  - z^)  Is  a 

characteristic  response  function  of  the  circuit  defined  In  the  follow- 
ing way.  With  reference  to  Figure  2-2  we  visualize  a hypothetical  gap 
cut  In  the  circuit  at  the  position  z^,  and  a voltage  of  one  volt 
applied  across  the  gap.  The  function  f(x  - z^)  Is  then  the  longitu- 
dinal component  of  the  electric  field  at  the  beam  position,  arising 
from  this  voltage.  In  the  one-dlmenslonal  theory  an  averaged  value  of 

f(x  - z ) across  the  beam  Is  used, 
o 

According  to  (2-5)  the  Induced  current  at  z^  Is  not  equal  to  the  beam 
current  at  z^  (possibly  modified  by  a constant  factor) , but  rather 
obtained  by  adding  weighed  contributions  from  the  beam  current  around 
the  position  z^.  Therefore,  the  Induced  current  Is  a smoothed,  or 
filtered  version  of  the  beam  current. 


The  Ramo  Equation  (2-5)  has  to  be  reformulated  to  fit  the  polarization 
model.  The  following  equation  applies: 


3I(Zo, 


^o> 


3z 


— ^ f f [x  +s(x,t) 
J " o ' o’  o' 


3s(x 


3t 


dx 


(2-6) 
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Figure  2-: 


*o 


X 


Schematic  representation  of  the  TWT  circuit 
and  beam,  and  the  definition  of  coupling 
function  f(x  - z ). 
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where  Is  the  undisplaced  coordinate  and  s(x^,  t^)  is  the  displacement, 
defined  by 


s(x^,  - X - \ (2-7) 

Equation  (2-6)  is  valid  with  overtaking  included.  The  set  of  indepen- 
dent variables  x^  and  t^  represent  the  undisplaced  coordinates,  l.e., 
the  reference  position  and  time  for  an  unmodulated  electron.  Note  that 
in  (2-6)  the  response  function  has  to  be  taken  at  the  displaced  posi- 
tion x^  + s(x^,  t^).  Hence,  the  equation  is  non-linear,  representing 
exactly  the  non-linear  beam-to-circult  coupling  discussed  earlier. 

The  integration  over  x^  takes  place  over  a sufficiently  large  interval 
around  to  include  all  non-zero  values  of  the  response  function. 

2.2  ELECTRONIC  EQUATIONS 

In  the  polarization  model  the  large  signal  electronic  equations  have 
the  form 


3u(z^, 


^o> 


at 


+ V 


9z 


- -nE[z„  + s(z.  t ), 


(2-8) 


as(z  , t ) 
o*  o 

at 


+ V 


az 


(2-9) 


The  electric  field  E in  the  driving  term  in  (2-8)  is  taken  at  the  dis- 
placed beam  position  rendering  this  term  non-linear.  The  field  consists 
of  the  circuit  field  E^  at  the  beam  position  and  the  space  charge 

field  E . 
s 
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E[Zo  + s(z^.t^),t^]  = E(,[2^  + s(Zo,t^),tJ  + E^[z^  + s(z„,t„),t„] 


o’  o' 


(2-10) 


These  are  given  by  the  following  expressions,  respectively: 


-I 


^'1^0  ^^*0**^0^ 


(2-11) 


nE  [z  + s(2  ,t  ),t  ] 
so  o o o 


^21  f R(k)^  ■*■  ' 

-J“p  u J -r- ' 

k=-«i 


oe 

I 


X 

o 


jk[x  + S(x  ,t  )1 

e dx  dk 

o 


(2-12) 


In  (2-11),  the  function  f[r  + s(z  t ) - x ] is  the  same  response  func- 

w 0 0 o 

tion  as  in  (2-6),  expressing  the  simple  fact  that  reciprocity  applies 
with  regard  to  coupling  between  beam  and  circuit. 


In  (2-12)  the  left  side  is  the  longitudinal  space  charge  field  averaged 

2 

over  the  beam  cross  section,  with  the  plasma  reduction  factor  R(k) 
specified  by: 


, , 21, (kb)  G,_(kb,  ka) 

R(k)  ■ R(kb,  ka)  - 1 ^ 


(2-13) 
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where 


G^Q(kb,  ka)  - I^(kb)  K^(ka)  + Kj^(kb)  I^(ka)  (2-14) 

This  Is  the  plasma  reduction  factor  for  a beam  constrained  to  motion  in 
the  longitudinal  direction  only  (confined  flow)  assuming  constant 
current  density  across  the  beam. 

In  the  subsequent  analysis  it  is  convenient  to  operate  with  the  Fourier 
transform  f(k)  of  the  response  function  f(z^).  By  definition 

f(Zo)  ~ 2^  I ° (2-15) 

k--» 


so  that 


f'(Xo  + S(x^,  t^)  - z^) 


00 

hf 


-jk[x  +s(x  ,t  )-z  ] 
kf(k)e  ° dk 


(2-16) 


The  Fourier  transform  f(k)  appearing  in  these  equations  is  obtained  by 
solving  the  appropriate  field  problem  for  the  configuration  shown  in 
Figure  1-2.  The  response  function  averaged  over  the  beam  cross  section 
has  the  Fourier  transform: 

21, (kb)  . 

f (k)  - f(kb.  ka)  . -ig—  (2-17) 

o 

The  preceding  discussion  shows  that  the  overall  set  of  circuit  and  beam 
equations  are  expressed  in  terns  of  two  characteristic  functions  of  the 
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wave  number  k,  namely  the  plasma  reduction  factor  R(k)  defined  in  (2-13), 
and  the  coupling  coefficient  f(k)  defined  by  (2-17).  In  contrast  to 
usual  TWT  analyses  the  present  theory  accounts  for  the  fact  that  the 
wave  number  k is  complex,  which  implies  that  R(k)  as  well  as  f (k)  are 
also  complex.  Hence,  there  are  small  but  noticeable  differences 
between  the  basic  models  used  in  the  present  theory  and  usual  TWT 
theory.  These  differences  show  up  even  in  the  small  signal  results, 
with  small  modifications  in  the  gain  as  well  as  synchronization 
conditions. 

2.3  INTRODUCTION  OF  NORMALIZED  VARIABLES  AND  EXPANSION  IN  HARMONICS 


The  following  normalized  variables  are  used  in  the  theory.  Independent 
normalized  variables: 


z = z 

o e o 


T = ojt 
o o 


6 ’ 

e 


where  u)  is  the  fundamental  frequency  and  where 


e V. 


The  normalized  rf  variables  of  the  circuit  are: 


v(z  , T ) 
o’  o 


^(-o’  ‘^o> 


(2-18) 

(2-19) 

(2-20) 


(2-21) 


(2-22) 
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o o 


I(z_,  t^) 


(2-23) 


and  those  of  the  beam: 


u(z  , t ) 

U<^o*  "o>  “ — 

o 


(2-24) 


S(Z„,  T ) • ^ sU,  t ) 
o o e o o 


(2-25) 


In  addition  to  the  normalized  displacement  S(Z  , T ) the  following 

o o 

variable  Is  Introduced: 


«(K.  2„.  I^) 


-iKSCZ^,  T^) 


(2-26) 


The  Inverse  relation  Is: 


-jKS(Z  T ) 

s ° ° - 1 - jK*(K,  2^  T^) 


(2-27) 


Apparently,  for  small  signals: 


Urn  ♦(K.  Z^.  T^)  - S(Z^.  T^) 
S-*o 


(2-28) 


The  function  ♦ (K,  Z^,  T^)  Is  seen  to  represent  the  non-linear  effects  in 
the  circuit  equation  as  well  as  the  electronic  equation. 


Since  the  rf  variables  are  periodic  in  time,  they  can  be  expanded  In 
harmonics  as  follows: 
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CO 


(2-29) 


V(Z. 


jiiT 


n o 


i(Z  , T ) 
o’  o 


i (Z  )e 
n o 


JnT 


(2-30) 


U(Z 


, T ) - 7 U 
o*  o n 


(Z^)e 


(2-31) 


■'o’  ■ S 


jnT 


(2-32) 


^o’  ” 2L,  ^2-33) 

—CO 

A A 

In  these  equations  v_ (Z  ) , i (Z  ),  etc  are  the  complex  amplitudes  of 

no  no 

the  nth  harmonics. 

Substitution  of  these  variables  into  the  circuit  equations  (2-3)  and 
(2-A)  and  the  electronic  equations  (2-8)  and  (2-9) , followed  by  separa- 
tlon  of  these  Into  harmonic  components,  results  In  the  overall  basic 
large  signal  equations  specified  In  the  next  section. 

2. A HARMONIC  SET  OF  LARGE  SIGNAL  EQUATIONS 

Following  the  procedure  outlined  In  the  preceding  discussion,  we  obtain 
the  following  set  of  large  signal  equations.  The  circuit  equations  for 
the  nth  harmonic  component  are  given  by: 
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Linear 


o 


(—\  jn»c  i - 0 


Linear 


--^+7r7 — ■ " I"  f ^ 

dZ  G Z-  \v_  I n 2ir  I 
o o Cn  \ Cn/  „ 


S (X  )e  dX  dK 
n o o 


1 

Non- 

linear 

1 

1 

q S— 00  00 

coupling 

1 

term 

1 

00 

r * A 

1 

1 § (X  )$  (-K,  X )e  ° dX  dK 

1 q o n-q  o o 

The  linear  part  of  the  coupling  term  is  retained  on  the  left-hand  side 
of  (2-35),  whereas  the  non-linear  right-hand  coupling  turn  is  identi- 
fied by  the  dotted  outline.  The  parameters  and  k appearing  in  the 

w n 

circuit  equations  are  defined  by: 
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Again,  the  non-linear  part  of  (2-39)  is  Identified  by  the  dotted  out- 
line. The  first  non-linear  term  arises  from  the  non-linearities  in  the 
beam  to  circuit  coupling.  The  second  and  third  terms  are  associated 
with  non-linearities  in  the  space  charge  field. 

In  these  equations  the  parameter  is  the  normalized  plasma  frequency 
of  the  beam 


b) 

_E 

u 


(2-41) 


The  non-linear  parts  of  the  circuit  equation  (2-35)  and  the  electronic 
equation  (2-39)  are  expressed  as  superpositions  of  cross-coupling 

A A 

terms  between  the  various  harmonics  of  the  rf  variables  v (z  ),  S (z  ), 
^ q'  o'’  q'‘  o" 

and  <>q(K,  z^).  The  dominant  cross-coupling  terms  are  those  centered 
around  the  fundamental  components  q = +1  and  q ° -1.  As  q deviates  more 
and  more  from  these  values,  the  cross-coupling  terms  become  smaller  and 
smaller  and  are  finally  negligible.  Expressed  differently,  it  suffices 
to  Include  only  a certain  number  of  harmonics.  The  exact  number  can 
be  determined  either  from  the  computational  procedure,  or  by  making  a 
reasonable  choice  based  on  previous  experience  from  corresponding 
klystron  programs  [1,2].  In  the  present  work  we  shall  take  q^^^  (or 
n ) equal  to  four.  The  fact  that  the  circuit  impedance  is 
vanishingly  small  when  the  harmonic  number  n exceeds  say  two  or  three, 

depending  on  the  helix  design,  implies  that  the  corresponding  circuit 
T,  in  \ 7 i7  \ 4 ru% 

^ i.  vci^xauxco  v / oiiVA  ^ \ *^  / • v/ka  s/ «.&»««.  itcAiAsi  wtt« 

n'  o n o e»  o 

A A A 

corresponding  dynamic  variables  ^o^ 

necessarily  negligible,  a fact  which  is  easily  observed  from  (2-39)  and 
(2-40).  The  higher  harmonics  in  the  beam  are  excited  predominantly  by 
the  fundamental  circuit  field,  and  to  a lesser  extent,  by  the  non- 
linear space  charge  field. 


It  is  noted  that  the  basic  large  signal  equations  are  dlf ferentlal- 
Integro  equations,  for  the  small  signal  as  well  as  the  large  signal 
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equations.  This  formulation  is  not  a consequence  of  the  polarization 
model  as  such,  but  is  rather  due  to  the  adopted  model  of  beam-to- 
circuit  coupling  mechanisms  discussed  in  Chapter  2.1.  The  dif ferential- 
integro  form  gives  rise  to  additional  beam  loading  terms  and  circuit 
loading  terms  that  do  not  appear  in  conventional  TWT  theory.  Mathemat- 
ically, the  difference  manifests  itself  in  the  appearance  of  the  plasma 

2 

reduction  factor  R(K)  and  the  coupling  coefficient  f(K),  which  in  the 
present  formulation  are  complex  rather  than  real.  It  should  be 
pointed  out,  however,  that  the  imaginary  parts  of  these  parameters  are 
quite  small  compared  to  the  real  parts,  so  that  the  corresponding 
corrections  are  quite  small.  However,  they  are  noticeable,  particu- 
larly as  far  as  optimum  synchronization  conditions  are  concerned. 

The  set  of  large  signal  Equations  (2-34)  - (2-40)  are  equally  valid 
for  small  signals,  which  is  obtained  as  a special  case  by  putting 
the  non-linear  terms  contained  within  the  dotted  outlines  equal  to 
zero.  The  major  steps  in  the  small  signal  case  will  first  be  treated 
separately,  as  it  forms  a necessary  part  of  the  subsequent  large  signal 
study. 


2.5  SMALL  SIGNAL  EQUATIONS 


According  to  the  foregoing  discussion  the  small  signal  dif ferential- 
integro  equations  are  specified  as  follows: 
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dS 

+ jn  S - U = t' 
dZ  ■'  n n 
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(2-45) 


As  noted  earlier  the  appearance  of  double  integrals  in  the  small  signal 
equations  is  a consequence  of  the  adapted  TVJT  model.  However,  due  to 
the  linearity  of  the  equations  they  can  be  solved  using  standard 
procedures. 

2.5.1  Small  Signal  Solutions  for  n 0. 

The  small  signal  solution  consists  of  a superposition  of  four  consti- 
tuent solutions.  Let  these  be  numbered  by  the  superscript  i,  where 
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i 1,  2,  3,  4.  A particular  rf  variable,  ay  S (z  ) is  then  given  by 

no  o y 

the  following  expression: 


S (Z  ) 
n o 


■L 


gd)  ^ n o 
n 


n 0 


(2-46) 


i=l 


.(i) 


where  is  the  normalized  propagation  factor  of  the  ith  solution. 

It  is  related  to  the  actual  propagation  factor  by: 


,(i) 


R(i)  _ 

- 6 


(2-47) 


The  amplitude  S is  independent  of  Z . If  we  consider  one  of  the 
n o 

independent  four  modes  specified  by  its  i“value  and  evaluate  the  cor- 
responding double  integrals  in  (2-43)  and  (2-44)  we  obtain; 
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By  this  procedure  we  find  that  the  1th  mode  of  the  small-signal  set 
(2-42)  - (2-45)  Is  specified  by  the  matrix  differential  equation: 


° 1 = 1,2, 3,4 


(2-51) 


where  the  column  vector  represents  the  rf  variables. 
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(2-52) 


The  4x4  system  matrix  Is  given  by 

n 
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n 0 

i = 1,  2,  3,  A 

Tlie  solution  of  the  small  signal  equation  (2-51)  is  straight  forward 
and  given  by: 


— i R Z 

S (Z  )(^>  = e ^ " °q 
~n  o 2n 


n ^ 0 (2-5A) 

i = 1,  2,  3,  A 


where  A is  the  normal  mode  amplitude  of  the  ith  mode,  and  where  q 
" n 

is  the  corresponding  state  vector.  The  overall  small  signal  solution 

is  a superposition  of  the  form  solutions  (2-5A) . 
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(2-55) 


where  the  state  matrix  Q 

^n 


consists  of  the  four  state  vectors. 
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n 0 


(2-56) 


q(2)  q(3) 
- n ~ n 
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Furthermore,  D is  a diagonal  matrix  containing  the  propagation  factors. 
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(2-57) 
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The  vector  A*  consists  of  the  four  mutually  independent  mode  amplitudes 

/ • \ ~ n 

( i ) 

, which  are  independent  of  Z^, 


A = 
- n 


Y 


(2-58) 


Equation  (2-55)  is  the  complete  and  general  small  signal  solution, 
expressed  as  a superposition  of  four  independent  normal  modes.  The 
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elements  of  and  are  obtained  from  the  appropriate  eigenvalue 
relation,  which  follows  from  (2-51).  For  the  ith  mode: 


n / 0 

i = 1,  2,  3,  4 


(2-59) 


The  propagation  factors  ^ are  obtained  by  putting  the  determinant 
equal  to  zero. 


l|  = 0,  n 0 (2-60) 

1 = 1,  2,  3,  4 

In  expanded  form  the  the  determinant  equation  is  given  by: 


n ^ 0 
i = 1,  2,  3,  4 

(2-61) 

The  four  solutions  of  this  equation  specify  the  four  propagation  factors 
of  the  nth  harmonic  component.  The  corresponding  state  vectors 
are  given  by: 
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n 0 

i = 1,  2,  3,  4 


(2-62) 

2.5.2  Small  Signal  Solutions  for  n = 0 

The  previous  equations  are  valid  when  n 0.  For  n = 0 the  system  is 
degenerate,  and  this  case  must  therefore  be  treated  separately.  Equa- 
tions (2-42)  - (2-45)  are  valid  for  n = 0 but  the  double  integrals 
(2-45)  - (2-50)  take  different  forms.  The  overall  solution  for  n = 0, 
corresponding  to  (2-55)  for  n 0,  is  given  by: 

So<^o>  ■ 8o(Zo>  i . <2-“> 

where  the  state  matrix  Q (Z  ) is  a function  of  distance  Z along  the 

~o  o o 

circuit. 
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(2-6A) 


The  normal  mode  vector  specifies  the  four  mutually  independent  time- 
average  or  "dc"  modes  corresponding  to  n = 0 


A 

-o 


(2-65) 


-'(1)  ■'(4) 

where  four  constants  A^  to  A^  represents  the  independent  normal 
node  amplitudes  of  the  n = 0 solution. 


In  normal  small  signal  operation  the  n = 0 mode  as  well  as  the  harmonic 
modes  corresponding  to  n are  not  excited.  However,  they  are 
excited  at  large  signal  levels  through  non-linear  effects.  Moreover, 
they  can  also  be  excited  at  small  signals  as  well  as  latge  signals 
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through  harmonic  injection  at  the  input  end  of  the  helix.  In  the 
present  treatment  we  shall  include  harmonic  injection  for  n = 0,  i.e., 
a "dc  injection  scheme".  The  motivation  for  this  is  that  we  like  to 
maintain  generality  and  not  exclude  the  possibility  of  improving  the 
Tl'JT  performance  by  "dc  injection",  which  would  supplement  the  usual 
second  harmonic  or  higher  harmonic  injection  schemes. 


For  this  purpose  it  is  appropriate  to  point  out  the  physical  signif- 
icance of  the  four  n = 0 modes,  each  of  which  is  specified  by  the  cor- 
responding column  in  the  state  matrix  9q(Z^)  in  (2-64). 


The  first  mode,  with  arbitrary  normal  mode  amplitude  determined 

by  the  excitation  conditions,  is  characterized  by  the  first  column  in 

(2-64).  Apparently,  this  mode  corresponds  to  a physical  situation  in 

“(1) 

which  the  circuit  carries  a small  dc  current  i “A  (the  second  ele- 

o o 

ment  in  the  column) . This  current  gives  rise  to  a dc  circuit  voltage 
proportional  to  the  circuit  resistance  and  to  distance  along  the 
circuit.  The  constant  voltage  gradient  arising  from  the  linear  variation 
of  voltage  along  the  circuit  serves  as  a driving  term  in  the  dynamic 
equations,  causing  a linearly  varying  dc  velocity  component  expressed 
by  the  third  element  in  the  column.  The  linearly  varying  velocity  com- 
ponent finally  gives  rise  to  a displacement  which  is  proportional  to 
the  square  of  distance  along  the  circuit,  as  shown  by  the  fourth  element 
in  the  column. 


In  summary,  this  first  dc  mode  has  a clear  physical  significance,  but 
is  hardly  of  much  practical  significance.  Its  effect  on  the  large 
signal  dynamics  would  come  from  the  constant  accelerating  (or  decelerat- 
ing) dc  field  arising  from  the  combination  of  Injected  dc  circuit  current 

and  the  circuit  dc  resistance  R . Aside  from  the  difficulties  in  inject- 

o 

ing  dc  current,  thermal  dissipation,  etc,  the  resistance  would  have  to 
be  fairly  high  in  order  for  this  effect  to  be  noticeable.  The  high  dc 
resistance  would  probably  be  detrimental  for  the  rf  operation.  Hence, 
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it  is  unlikely  that  a dc  injection  scheme  based  on  the  first  mode  is 

of  any  practical  value.  Moreover,  the  excitation  of  this  mode  directly 

from  the  beam  due  to  non-linear  effects  under  large  signal  conditions 

is  also  prohibited,  provided  special  arrangements  are  not  made  to 

allow  a dc  current  to  flow  in  the  circuit.  This  is  usually  not  the  case. 

Therefore,  in  the  large-signal  calculations  we  shall  disregard  this 

''  (1) 

mode  altogether,  putting  = 0. 


The  second  dc  mode,  corresponding  to  the  second  column  of  the  state 
matrix  current  in  the  circuit  but  a constant  dc  voltage 

A ^ ° (2)  , 

V =A  . The  voltage  corresponds  to  a dc  velocity  change  of  A /2 
o o . o 

and  a linearly  changing  displacement  A^  Z^/2,  as  seen  from  the  two 
last  elements  in  the  column.  This  mode  could  conceivably  be  used  in  a 
"dc  injection  scheme".  For  instance,  it  would  characterize  a physical 
situation  in  which  the  two  circuit  sections  of  a normal  TWT  are  not  at 
the  same  dc  voltage.  The  difference  in  voltage  corresponds  to  the 
normal  mode  amplitude  of  the  second  dc  mode.  As  shown  in 

Section  2.7.2  the  second  dc  mode  is  not  excited  by  non-linear  effects 
in  the  beam-to-circuit  coupling,  but  only  from  "dc  injection"  directly 
on  the  circuit. 


The  third  and  fourth  dc  modes  are  characterized  by  zero  circuit  voltage 
and  current,  as  seen  directly  from  (2-64).  The  third  dc  mode  specifies 

A A (3) 

a constant  change  in  relative  dc  velocity  given  by  , asso- 

ciated with  a linearly  increasing  or  decreasing  displacement.  Phys- 
ically, this  mode  can  be  visualized  as  arising  from  a small  change  in 
dc  velocity  v^.  It  follows  from  the  subsequent  large  signal  theory 
that  the  mode  basically  describes  the  process  of  slowing  down  of  the 
electrons  under  large  signal  operation. 

The  fourth  dc  mode  represents  a constant  displacement  with  no  change  in 
velocity.  The  physical  significance  of  this  mode  is  less  obvious.  It 
can  conceivably  be  excited  through  some  kind  of  time  delays  arising 


26 


from  non-linear  effects  under  large-signal  operation.  Regardless  of 
physical  interpretations,  the  third  and  fourth  dc  modes  corresponding 
to  n = 0 are  certainly  excited  under  large-signal  conditions,  as  shown 
in  Section  2.7.2,  and  must  be  included  as  part  of  the  overall  mode 
configuration  in  the  large-signal  analysis. 

2.6  COMPARISON  WITH  CONVENTIONAL  TWT  THEORY 


In  the  conventional  TWT  model  the  circuit  is  essentially  visualized  as 
being  at  the  same  position  as  the  electron  beam,  with  an  appropriate 
reduction  in  circuit  impedance  depending  on  the  frequency.  The  fact 
that  the  reduction  depends  on  frequency  rather  than  wave  number  gives 
rise  to  certain  deficiencies  and  shortcomings  in  the  conventional  TWT 
model.  We  shall  briefly  discuss  these  and  point  out  how  these  short- 
comings are  eliminated  in  the  present  model. 

The  differences  between  the  two  models,  at  least  at  small  signal  level, 
have  nothing  to  do  with  the  use  of  polarization  variables,  but  arise 
from  a different  description  of  the  coupling  between  beam  and  circuit. 
In  the  present  model  the  coupling  is  essentially  specified  by  the 
coupling  factor  f which  appears  in  the  off-diagonal  coupling 

terms  in  the  overall  system  matrix  M in  (2-53).  The  essential  point 
is  that  the  argument  in  the  f ” function  is  the  actual  propaga- 
tion factor  of  the  ith  mode.  In  conventional  theory  the  argument  is 
taken  to  be  equal  to  or  for  the  nth  harmonic,  Hence 

the  coupling  factor  is  taken  to  be  the  same  for  all  four  modes,  and 
is  frequency  dependent,  rather  than  wavelength  dependent  as  in  the 
present  model.  At  first  glance  this  difference  may  seem  fairly  insig- 
nificant, because  the  four  propagation  factors  do  not  deviate  much  from 
the  cold  circuit  propagation  factor  6^.  This  is  true  under  normal 
operation  where  the  present  model  yields  only  slightly  different  gain 
response,  although  the  optimum  synchronization  condition  is  signifi- 
cantly different.  However,  the  differences  first  come  to  proper  light 
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when  circuits  with  heavy  loss,  i.e.,  attenuators,  are  considered.  In 
such  cases  there  exists  an  essential  disagreement  between  conventional 
TWT  theory  and  the  present  theory. 

As  the  circuit  loss  is  increased,  the  loss  parameter  defined  in 

(2-37)  becomes  larger  and  larger.  From  consideration  of  the  wave 

equations  (2-59)  - (2-62)  it  is  easy  to  show  that  in  the  limit  of  very 

large  k , the  circuit-to-beam  coupling  disappears  for  the  two  circuit 
n 

modes.  These  are  specified  by  the  propagation  factors 


lim 


2 


(2-66) 


The  coupling  factor  f 
ment  (see  (2-17)). 


becomes  negligible  due  to  the  large  argu- 


lim 

n 


= 0 


(2-67) 


Hence,  it  follows  from  (2-53)  that  the  off-diagonal  coupling  terms 
disappear,  rendering  the  system  completely  decoupled  from  the  beam  as 
far  as  the  heavily  damped  circuit  modes  are  concerned.  This  is  not 
predicted  by  the  conventional  TWT  model  in  which  the  coupling  factors 
are  taken  to  be  Independent  of  the  propagation  factors. 


The  remaining  two  modes  in  a lossy  circuit  section  are  essentially  Che 
two  space  charge  waves  in  the  beam,  which  have  a finite  non-zero  coup- 
ling to  the  lossy  circuit  because  the  propagation  factors  are  finite, 
and  therefore  f different  from  zero. 

In  conclusion,  the  present  model  distinguishes  itself  by  its  ability  to 
account  properly  for  arbitrarily  large  loss  in  the  circuit  structure 
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This  is  a great  convenience  in  numerical  calculations  of  multi-section 
HTT's,  where  the  same  basic  equations  can  be  used  regardless  of  their 
loss. 

In  addition  to  this  discussion  we  shall  be  a little  more  specific  about 
the  fundamental  difference  between  the  conventional  Pierce  model  and 
the  present  model.  It  is  convenient  to  use  the  dispersion  relation 
(2-61)  as  a basis,  since  it  is  easily  converted  to  the  Pierce  form 
by  changing  the  coupling  factor  f and  the  plasma  reduction 

factor  r(b<«)  to  the  approximations  used  in  the  Pierce  model.  In  the 
Pierce  model  the  small  signal  TOT  performance  is  completely  described 
by  specifying  four  independent  parameters:  the  gain  parameter  C,  the 
loss  parameter  d,  the  velocity  parameter  b,  and  the  space  change  param- 
eter QC.  Once  these  are  specified,  the  gain  characteristics  are  com- 
pletely determined  regardless  of  other  details  of  circuit  and  beam  param- 
eters. 


The  four  basic  Pierce  parameters  are  related  to  the  parameters  used  in 
the  present  theory  by  the  following  equations: 


C = 


^ G Z-f 
4 O C 


1/3 


d 


1 

C 6 


e 


(2-68) 
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where  we  consider  only  the  fundamental  frequency,  In  (2-68)  a is  the 
loss  factor  of  the  helix  (imaginary  part  of  the  complex  cold  circuit 
propagation  factor).  It  is  related  to  the  loss  factor  k defined  in 
(2-37)  by  the  equation: 


K = 


(2-69) 


The  argument  v /v  in  the  coupling  coefficient  f (v  /v  ) indicates  that 

O ^ O L« 

in  the  Pierce  theory  the  normalized  propagation  factor  B = B/B^  is  put 
equal  to  Bp/B  = v /v  . This  means  that  the  propagation  factor  is  put 

0 6 O O 

equal  to  the  cold  circuit  propagation  factor  B_,  thus  ensuring  that  the 

o 

gain  parameter  C in  (2-68)  is  evaluated  on  the  basis  of  the  cold  inter- 

2 

action  impedance  Z f (v  /v  ) at  the  beam  position. 

O 0 0 

The  argument  of  unity  in  the  plasma  reduction  factor  R(l),  appearing  in 
(2-68),  indicates  that  in  the  Pierce  theory  R is  evaluated  by  putting 
B = 6/Bg  = 1,  i.e.,  B = B^. 


Inversion  of  the  equations  specified  by  (2-68)  yields  the  following 
relations: 


K = 


o C 


1/2 


= X;^(Cd,Cb) 


V 

= 1 + Cb  = X2(Cb) 


(2-70) 


2 2 
n R(l)  = 
P 


2 = X3(C,QC) 
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Let  us  now  express  the  dispersion  relation  (2-61)  in  terms  of  the  four 
Pierce  parameters  by  substituting  from  (2-69)  and  (2-70)  into  (2-61), 
putting  n = 1 and  for  simplicity  omitting  the  superscript  i.  The 
result  is 


B 


4 


3 2 

+ B 


1 - X3(C,QC) 


R(B)^ 

R(l)^ 


(2-71) 


X2(Cb)‘ 


Xj^(Cd,Cb) 


1 - 


X3(C,QC) 


R(B)^ 

R(l)^ 


= 0 


This  is  the  small  signal  dispersion  relation  valid  for  the  present 
model,  but  expressed  in  terms  of  the  four  Pierce  parameters  defined  in 
(2-68).  It  is  noticed,  however,  that  in  addition  to  these  parameters, 
two  more  parameters  are  needed,  namely: 


a)  Modified  coupling  parameter: 


4 


1 


f(B)^ 

i(x2^Ch)y 


b)  Modified  space  charge  parameter: 


4 


2 


R(B)^ 

R(l)^ 


(2-72) 


(2-73) 
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The  essential  difference  between  the  present  model  and  Pierce's  model 
can  be  appreciated  from  this  exposition.  In  the  Pierce  model  the 
factors  and  are  both  approximated  by  putting  them  equal  to  unity. 
Inspection  of  the  dispersion  relation  (2'-71)  in  this  approximation 
reveals  that  the  four  propagation  factors  are  specified  fully  by  the 
four  Pierce  parameters  C,  d,  b,  and  QC.  As  an  example,  this  model 
predicts  identical  gain  response  for  two  TWT  designs  with  different 
circuit  radii  (and  possibly  different  beam  radii  as  well)  provided  the 
four  basic  Pierce  parameters  are  the  same  in  both  cases.  And  this  may 
very  well  be  possible,  since  a larger  circuit  diameter  can  be  offset  by 
a basically  larger  interaction  impedance  measured  at  the  circuit,  so  that 
the  interaction  impedance  is  the  same  at  the  position  of  the  beam,  i.e., 
the  gain  parameter  C is  the  same. 

In  the  present  model  the  four  basic  Pierce  parameters  must  be  supple- 
mented by  the  parameters  and  defined  in  (2-72)  and  (2~73’>.  These 
are,  by  themselves,  functions  of  B,  or  more  specifically,  functions  of 
6b  and  6a,  which  can  be  expressed  as  B(6^b)  and  B(6^a) , respectively. 
Hence,  as  far  as  dependence  on  TV7T  parameters  are  concerned,  and  ^2 

are  functions  of  the  normalized  beam  and  circuit  radii  6 b and  6 a. 

e e 

Hence,  we  reach  the  conclusion  that  our  TWT  model  requires  six  basic 

parameters  rather  than  the  four  required  in  the  Pierce  model.  The 

additional  parameters  are  the  normalized  beam  and  circuit  radii  6 b and 

e 

6^a.  These  are  not  requ'*’ — in  Pierce's  model. 

As  an  example  consider  two  TWT  designs  both  having  the  same  Pierce 
parameters  C,  d,  b,  and  QC,  but  different  normalized  beam  and  circuit 
radii.  The  Pierce  model  gives  identical  gain  characteristics  for  these 
two  designs,  whereas  the  present  model  predicts  a difference  in  gain 
depending  on  the  normalized  beam  and  circuit  radii  in  the  two  designs. 
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The  circumstance  that  the  two  factors  and  ^ particular  mode 

depend  on  the  actual  propagation  factor  g as  well,  leads  us  to  the  con- 
clusion that  the  differences  in  gain  characteristics  obtained  from  the 
two  models  will  be  larger,  the  more  the  complex  propagation  factor  B 

deviates  from  B and  B„.  Therefore,  in  high-gain  tubes  and  in  attenu- 
e C 

ating  sections  we  expect  fairly  large  differences  between  the  two 
models. 

2.7  LARGE  SIGNAL  SOLUTION 


The  large  signal  equations  are  given  by  (2-34)  - (2-39)  which  also 
contained  the  small  signal  equations  discussed  as  a special  case  in  the 
preceding  sections.  The  motivation  for  solving  the  small  signal  case 
in  detail  is  not  dictated  merely  from  Interest  in  small  signal  response, 
but  more  from  the  fact  that  certain  small  signal  parameters  are  needed 
in  the  subsequent  large  signal  solutions. 

In  view  of  the  fact  that  the  left  hand  sides  of  the  general  large  signal 
equations  are  Identical  with  the  small  signal  equations  (2-42)  - (2-45), 
we  can  take  the  general  large  signal  solution  to  be  essentially  of  the 
same  form  as  the  small  signal  solution  (2-55)  and  (2-63),  namely: 


a (Z  ) 

= Qe 

° A (Z^) 

n 0 

(2-74) 

-n  o 

n 0 

'•O  0 

= Q (Z  ) A (Z  ) 

^o  o ~o  0 

n ==  0 

(2-75) 

The  state  matrices  and  QqCZ^j)  and  the  propagation  factors  are 
identical  in  the  small  and  large  signal  case.  The  large-signal  expres- 
sions (2-74)  and  (2-75)  differ  from  the  small  signal  expressions  only 

A A 

In  the  fact  that  the  normal  mode  amplitude  vectors  A^(Z^)  and 
are  taken  to  be  functions  of  distance  In  the  small  signal  case 

these  vectors  are  constant.  The  variations  of  the  normal  mode 
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amplitudes  are  caused  by  the  non-linear  interactions  between  modes. 

These  manifest  themselves  as  saturation  effects  in  the  fundamental 
component  (the  corresponding  vector  A^(Z^)  becomes  smaller  as  the 
signal  amplitude  grows  along  the  circuit),  and  in  non-linear  excitation 
of  higher  harmonics  and  intermodulation  components  (the  corresponding 
vectors  A^(Z^)  grow  along  the  tube). 

The  essence  of  the  large  signal  calculations  is  therefore  to  determine 
the  variations  of  the  normal  mode  vectors  A (Z  ) and  A (Z  ) along  the 
length  of  the  tube  for  the  fundamental  component  and  all  harmonics  of 
interest.  Once  this  has  been  achieved  the  general  large  signal  solution 
for  these  nth  harmonic  is  given  by  (2-74) , and  for  the  time  average 
component  by  (2-75).  The  overall  large  signal  solution  is  obtained  by 
superposing  the  harmonics: 


^ ° + y a (Z„)e'  ° 


n =— <xj 


'O  o o o 


'O  o /•  ^ ~n ' o' 
n ^0 


(2-76) 


E-jD  Z ^ 

Q e ° A 
~n  -n 


-jD_Z„  ^ jnT 

(Z^)e  ° 

o 


n^O 


It  should  be  emphasized  that  the  assumed  form  of  the  large  signal 
solutions  (2-74)  and  (2-75)  are  not  approximations,  but  perfectly 
valid  at  all  signal  levels.  The  normal  mode  amplitudes  A^(Z^)  and 
A^(Z^)  contain  the  necessary  large  signal  information  of  amplitude  as 
well  as  phase  variations  of  all  four  modes  and  all  harmonics. 
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2.7.1  Differential  Equations  for  the  Large  Signal  Normal  Mode 
Amplitudes 

Mathematical  processing  of  the  large  signal  equations  (2-34)  - (2-39) 
results  In  the  following  differential  equations  for  the  normal  mode 
amplitudes  defined  in  (2-74)  and  (2-75): 


dA  (Z  ) jD  Z 
~n  o ~no  1 

= e Q F (Z  ) 

~n  ~n  ° 


dZ 


n 0 


(2-77) 


dA  (Z  ) , 

° = Q (Z  )"^  F (z  ) 

dZ  ~o  o -o  o 


n = 0 


(2-78) 


The  vectors  F (z  ) and  F (z  ) are  non-linear  driving  terms  which  are 
-n  o ~o  o 

zero  in  the  small  signal  case.  In  general: 


n o 


“I 


«2a("o> 


»3n^'o> 


n it  0 


(2-79) 


-o  o 


3o  o 


a “ 0 


(2-80) 
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and  appearing  in  these  vectors  will  be 


The  elements  H2 
discussed  later 


nl 


The  differential  equations  (2-77)  and  (2-78)  can  be  Integrated  directly, 
yielding: 


jD  X 
~n  o 


F (x  ) dx  + A (o) 
**’  Q o o - H 


n 0 


Z 

o 

Q (x  ) ^ F (x  ) dx  + A (o) 
-Jo  O -0  0 0 -0 

o 


n = 0 


(2-81) 


(2-82) 


In  these  equations  A (o)  and  A (o)  are  the  initial  values  of  the  normal 
-n  -o 

mode  vectors  at  = 0,  i.e.,  of  the  input  end  of  the  circuit.  Because 
of  the  low  signal  level  at  the  input,  these  are  identical  to  the  small 

signal  values.  Without  harmonic  Injection  A (o)  = 0 for  n ^ ±1.  With 

- a 

harmonic  Injection  the  corresponding  vectors  A^(o)  are  non-zero. 

2.7.2  Excitation  of  the  Time  Average  Modes  (g  = 0) 

The  0=0  mode  amplitudes  are  given  by  (2-82).  Evaluation  of  ^ 

and  use  of  (2-80)  yield  the  following  expressions  for  the  excitations 
of  the  "dc  modes": 


A^^\z  ) - A^^\o)  - constant 
0 0 o 


(2-83) 


A^^\z  ) » A^^^(o)  - A^^^  « constant 
0 0 o o 


(2-8A) 
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(2-85) 


r 

= -J 


Z 


^(4) 


r 


(2-86) 


According  to  these  equations  the  two  first  modes  do  not  depend  on  z^. 
They  are  not  excited  by  the  non-linearities  in  the  beam,  but  arise  from 
a possible  external  application  of  dc  current  or  dc  voltage  in  the 
circuit,  as  was  discussed  extensively  In  2.5.2.  On  the  other  hand,  the 
two  last  dc  modes  are  excited  by  the  non-linear  interactions  in  the 
beam.  As  described  in  2.5.2  they  describe  the  general  slowing  down 
process  taking  place  as  the  signal  amplitude  (and  therefore  ^^^(Z^)) 
increases. 


As  discussed  in  2.5.2  the  first  mode  is  zero  unless  the  circuit  is  excited 
with  a dc  current,  which  normally  is  not  practical.  Hence,  mode  2 is 
the  only  one  available  for  a possible  "dc  Injection"  scheme,  in  which 
an  additional  dc  voltage  equal  to  A^  (o)  is  applied  to  the  circuit. 

Using  (2-64)  and  (2-75)  in  conjunction  with  (2-83)  - (2-86)  the  actual 
time-average,  or  dc,  variables  are  obtained  as  follows: 

In  the  circuit  ("dc  injection") : 

V (Z  ) = v (o)  * a'  “ constant  (2-87) 

o o o o 


i (z„)  - 0 

o o 


(2-88) 
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In  the  beam: 


U (Z^)  = 1/2 
o o o 


- 


H_  (x  ) dx 
3o  o o 


(2-89) 


\ r ] r 

2 - J + J 


Z„  + / x^H3^(x^)  dx^  (2-90) 


The  relative  reduction  of  time-average  velocity  is  specified  by  U (z  ) 

o o 

in  (2-89),  where  the  first  term  is  due  to  the  extra  dc  voltage  of 

- (2) 

Aq  applied  on  the  circuit  ("dc  injection") , and  the  second  term  repre- 
sents the  slowing  down  process  caused  by  non-linear  effects  during 
the  signal  growth  along  the  beam. 


2.7.3  The  Non-Linear  Driving  Terms  in  the  Large  Signal  Equations 

The  non-linear  driving  terms  F^(Zq)  and  F^(Zq)  in  the  large  signal 

equations  (2-77)  - (2-82)  contain  the  elements  H,  (z  ) , H,  (z  ) , and 

zq  o oq  o 

H3^(z^),  which  are  obtained  from  the  basic  large  signal  equations 
(2-35)  and  (2-39).  The  summation  over  q is  eliminated  by  using  the 
following  relations: 


JK  ^ . n[5_,(-K,x^)  - S„(x^)]  (2-91) 


q —a 


00  A 


2^  ■ dT  L*n<-’'>’'o)  - S. («<,)]  (2-«) 

o o ^ ^ 

q s^oo 


Hereby,  the  expression  for  H-_(Z^)  simplifies  to 

Zb  O 
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00 

K=-“ 


Kf(K)e 


-jKZ^ 


(2-93) 


oo 

/ - S.ho)] 


jKx 

e ° dx  dK 
o 


X =-“ 

o 


which  is  valid  for  all  n,  including  zero.  Hence,  for  n = 0 


“2<,<2o>  - 0 


(2-94) 


which  we  have  already  made  use  of  in  (2-80) . 

Proceeding  to  H-  (z  ) we  find,  after  integration  by  parts: 


H-  iz) 
3n  o 


00  OO 

p=-oo 


-jKZ 


00 

/ 


v„(x  )e 
P o 


jKx^ 


dx  dK 


00 

X =-* 

o 


jKx 

e ° dx  dK 
o 


p«-00  I^-oo 


..  jKx 

((i  (-K,x  )e  ° dx  dK 
p o o 

(2-95) 


2,7.4  Transformation  of  the  Double  Integrals 


The  non-linear  coupling  terms  H2^  and  are  expressed  as  double 
Integrals  over  K and  x^.  These  are  not  easily  amendable  to  numerical 
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computer  calculations.  Therefore,  we  have  transformed  the  double 
integrals  to  a different  form,  which  is  derived  as  follows: 

Let  us  consider  a general  double  integral  of  the  form: 


00 

"'V-k  / 

K“-“ 


-jKZ 

a(K)e  °b(K 


00 

/ 


jKx 

C(K,x  )e  ° dx  dK 
o o 


X =-“> 

o 


(2-96) 


All  the  double  integrals  in  (2-93)  and  (2-95)  are  represented  by  this 
general  expression  provided  we  make  the  appropriate  identifications  of 

A A 

the  functions  a(K),  b(K,Z^)  and  C(K,x^).  For  instance,  the  function 

a(K)  represents  Kf(K)  in  (2-93),  Kf(K),  R(K)^,  or  KR(K)^  in  (2-95).  The 

function  6(K,Zp)  is  unity  in  (2-93),  it  is  unity  or  J _ (K,Z  ) in  (2-95). 

The  function  C(K,x^)  represents  several  different  variables,  such  as 

Tj  (-K,x  ) - S (x  )1,  V (x  ) , and  $ (-K,x  ). 

L n o n o J’  p o ’ p ’ o 

Regardless  of  which  variables  are  Involved,  the  function  a(K)  is  a 

A A 

slowly  varying  function  of  K.  The  functions  b(K,Z^)  and  C(K,x^)  are 
"wavelike"  functions,  almost  periodic  in  Z^  and  x^. 


Let  us  make  the  following  expansion  of  C(K,x^): 


-J ♦£<*„) 


(2-97) 


where  the  phase  angle  ^^(Xq)  is  chosen  so  that  the  preceding  factors 

A 

C^(Vi^,x^)  are  slowly  varying  functions  of  x^.  With  the  rf  variables 
represented  by  the  C(K,x  ) function,  this  Is  always  possible.  Expanding 

A 

the  slowly  varying  functions  C.(K  ,x  ) and  ^.(x  ) in  Taylor  series,  we 
can  convert  the  double  Integral  (2-96)  to  the  following  form: 
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T(z^) 


00 


r=0 


— [ 
dK’^  L 


a(K)b(K,2^) 


K - 


dz 


(2-98) 

This  is  an  exact  transformation  of  the  double  integral.  However,  since 
we  have  taken  care  to  expand  the  "wavelike"  function  C(K,x^)  in  the 
sum  (2-97)  80  that  Cj^(K,x^)  are  almost  independent  of  Z^,  the  deriva- 
tives of  this  function  are  negligible  and  can  be  disregarded  compared 
to  the  first  term  in  the  expansion  (2-98).  This  means  that  the  only 
term  to  be  retained  is  the  r = 0 term.  Hence 

(2-99) 

I 


where 

= (2-100) 
o 

The  approximation  (2-99)  of  the  exact  expression  (2-98)  corresponds  to 
disregarding  the  growth  or  decay  of  the  amplitude  of  the  function 
C(K,x^)  over  some  typical  "interaction  length".  The  interaction  length 
corresponds  to  the  nominal  width  of  the  coupling  function  between 

circuit  and  beam  (see  Figure  2-2),  if  a(K)  represents  £(K).  In  the  case 
that  a(K)  stands  for  R(K)2,  the  Interaction  length  is  typically  equal 
to  the  width  of  the  Green  function  for  the  space  charge  field. 
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Regardless  of  details,  (2-100)  is  a sufficiently  good  approximation 
for  evaluation  of  the  non-linear  coupling  integrals  (2-93)  and  (2-95). 
The  practical  problem  is  then  reduced  to  determination  of  the  functions 
C (K. ,x  ) and  the  phase  angles  (>.(x  ) in  the  expansion  (2-97).  This 

36  O J6  O 

information  is  extracted  from  the  appropriate  rf  variable  representing 
the  C(K,x^)  function,  as  discussed  earlier.  In  most  cases  the  ^q(K,x^) 
functions  are  involved  in  this  procedure,  and  a major  part  of  the 
remaining  details  is  concerned  with  the  expansions  of  the  various 
harmonics  of  the  <Ji(k,x^)  - function  defined  by  (2-26).  These  expansions 
are  necessary  details  for  the  numerical  procedure,  but  do  not  provide 
more  insight  into  the  architectural  structure  of  the  large  signal  theory 
and  program.  They  are  therefore  not  included  in  the  present  report. 

Nor  are  details  concerning  multi-section  TVJT's,  which  are  easily 
handled  by  the  present  theory,  and  considerations  of  launching  condi- 
tions and  extraction  of  the  rf  signal  at  the  output. 
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